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Abstract 


Identifiability of linear, time- invariant systems in physical coordinates is discussed 
in this paper. It is shown that identification of the system matrix in physical coordinates 
can be accomplished by determining a transformation matrix that relates the physical 
locations of actuators and sensors to the test-data-derived input and output matrices. 
For systems with symmetric matrices, the solution of a constrained optimization problem 
is used to characterize all the possible solutions of the transformation matrix. Conditions 
for the existence of a unique transformation matrix are established easily from the explicit 
form of the solutions. For systems with limited inputs and outputs, the question about 
which part of the system can be uniquely identified is also answered. A simple mass- 
spring system is used to verify the conclusions of this study. 


Introduction 


Modal updating has gained some research interest largely because the finite-element 

models of structural systems are usually not accurate enough for dynamic simulation and 

control design purposes. Most modal updating techniques were developed with the goal 

of using experimental data to improve the accuracy of finite-element models. Attempts 

include using constrained parameter optimization, nonlinear dynamic programming, and 

complicated iterative schemes to update the finite-element model so that it will fit a 
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given large set of test data. Similar concepts also have been used to develop damage 
detection techniques. The goal aims for “updating” a mathematical model of the un- 
damaged structure by using experimental data of the damaged structure and then, using 
the results to detect the locations of damaged structure elements. However, uniqueness 
of the solution has not been seriously investigated. If the uniqueness question is not 
answered, one will never be able to settle the dispute about whether model updating and 
damage detection techniques will really work in practical situations or not. 

This paper discusses identifiability of linear, time invariant systems with a limited 
number of inputs and outputs. The main conclusion is that if there are not enough 
actuators and sensors, the system model in physical coordinates can not be uniquely 
determined. Consequently, all existing damage detection techniques, no matter how 
sophisticated, can never exactly predict the locations of damaged parts of a structure. It 
is also pointed out in this paper that by knowing the physical locations of the actuators 
and sensors, it is possible to exactly identify a part of the physical- coordinates model. 
Therefore, if one s purpose is to update a physical model of the structure, it is important 
to place actuators and sensors at the right locations. 

The organization of this paper is as follows. First, the identification problem con- 
sidered here is described by a statement of the problem. It shows that the identification of 
system matrix can be accomplished by determining a transformation matrix. Then, iden- 
tifiability of systems with unsymmetric matrices is discussed. After that, a constrained 
optimization problem is solved with the results to be used in the discussion of identifi- 
ability of systems with symmetric matrices. The conclusions of this study are basically 
the same as those of Ref. [1], However, here we explicitly solve for all the possible solu- 
tions of the transformation matrix. Therefore, conditions for the existence of a unique 
transformation matrix can be established very easily from the form of the solution. The 
question about which part of the system is exactly identified is also answered. A simple 
numerical example is included. 
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Statement of Problem 


The system identification problem considered here is described as follows: Given 
an input-output transfer function matrix G(s) € R mxl (derived from test data) and an 
input matrix B € R nxl and an output matrix C € R mXn (derived from physical locations 
of actuators and sensors), determine a matrix A € R nXn such that 

C(sl- A)~'B = G{s) (1) 

We are interested in the uniqueness of the solution. If A is uniquely determined, then it 
represents the system matrix in physical coordinates. For model updating and structural 
redesign purposes, it is important to assure that the system matrix identified is indeed 
in physical coordinates. 

For a given transfer function, there is no unique state-space representation. Ex 
isting system realization techniques can find at least one state-space representation 
(A,B,C) that also satisfies G(s) = C(sl — A)~ l B. Since the two triples (A,B,C) 
and (A, B, C) describe the same system, they are related by a similarity transformation 

A = P~ X AP 

B = P~ l B ( 2 ) 

C = CP 

P € R nXn is called the transformation matrix. Therefore, identification of the A matrix 
in physical coordinates can be accomplished by determining the transformation matrix 
P, which relates the physical-coordinates model (A,B,C) to a given test-data-denved 
model (A, B, C). 

In the following sections, we will show that the transformation matrix is generally 
not unique if there is not enough sensors and actuators. Identification of first-order 
systems with unsymmetric system matrices will be discussed first. Then, a constrained 
optimization problem will be solved along with uniqueness of the solution studied. The 
results are used to discuss identifiability of symmetric system matrices and mass and 
stiffness matrices of structural dynamics systems. 
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Identification of Unsymmetric System matrix 


For general linear systems with unsymmetric system matrices, the transformation 
matrix can be determined from solving Eqs. (2b) and (2c), which are rewritten as 

B = PB 
C = CP 

Let the singular value decompositions (SVD) of B and C be 

B = , C = Uc'BcVq (4) 

Then, by premultiplying and postmultiplying Eqs. (3) by appropriate matrices and by us- 
ing the orthogonality property of singular value decompositions, Eq. (3) can be rewritten 
as 

Vq BVg = VJPUq Eg 
UqCUq = Y>cVq PUg 

Define a new matrix P = Vq PUg, then the above equations become 

VJBV b = PE B 
UqCUq = EcP 

Therefore, determination of P can be accomplished from determining P. 

Let P matrix be partitioned into 


(5) 


( 6 ) 


-[*> 

l Pn 


P\2 

Pi 1 


P n € R mxl 

Pn € R ( n - m )* 1 


Pn 6 R mx ^ n -‘) 

Pn e 


(7) 


Since B G R nxl and is a full-rank matrix, B has / nonzero singular values. That is, 
— 0 ’ w ^ ere [ ] is an / x / diagonal matrix with nonzero singular values 

on the diagonal. Therefore, Eq. (6a) uniquely determines the first / columns of P matrix. 
That is, 


Pn 
Pi i 


= vj BV s l(<T,) e 




(8) 


Similarly, because C matrix has m nonzero singular values, Eq. (6b) uniquely determines 
the first m rows of P matrix. That is 


Pu Pii}=[{a i ) c }~ l UlCU t 


(9) 


4 


There is only one partition, that is P 22 , that does not appear in the solutions, Eq. (8) 
and Eq. (9). If / > n and rank(B) = n, which implies the case that the number of 
actuators is greater than or equal to the number of states, then Eg is a full rank matrix, 
and hence, Eq. (6a) uniquely determines P. Similarly, if m > n and rank(C) = n, then 
Eq. (6b) uniquely determines P. If both / and m are less than n, then the solution of P 

is non-unique, because P 22 is undetermined. 

Therefore, the necessary and sufficient condition for the transformation matrix P 
to be uniquely determined is that either rank (B)=n or rank (C)=n. In other words, there 
must be at least as many sensors or as many actuators as the order of the system so that 
the A matrix can be uniquely identified. Also note that Ai is determined twice in Eqs. (8) 
and (9). The Ai determined from Eqs. (8) and (9) can be expressed, respectively, by 

Pu = f / J 1 flV / B[(<T,) S ] _1 , Ai = [(^i) c r lu cCU B i ( 10 ) 

where V& is the first m rows of Vj and U Bl is the first / columns of U B . The above two 
solutions of P n must be equal for consistency. In fact, 

= [( <T i)cY'VcCUB i — 4 [McWciBVb = UlCU Bx [{(Ti) B } 

— ► U c [{<n)c\ v ciBV B Vl = U c UZCU Bl [(<ri) B ] v 8 — > CB = CB 
The condition that CB = CB must be satisfied if (A,B,C) and (A,B,C) indeed repre- 
sent the same system. Therefore, if the two solution of P n in Eq. (10) are not consistent, 
then the given test-data-derived model (A, B, C) is a wrong model. 

Before we move on to discuss identification of systems with symmetric matrix, let us 
study the solution of a constrained optimization problem. Non-uniqueness of the solution 
of this optimization problem will prove to be useful in the discussion of identifiability of 

symmetric system matrices. 

A Constrained Optimization Problem 

The constrained optimization problem considered is: Given two matrices Q and 
R, both G R pxn and of ranks < n (p can be > or < n), and a positive-definite symmetric 

matrix M € R nXn , find all P € R nXn that 
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( 11 ) 


minimize J =|| Q - RP ||2 
subject to P t MP = I n 
This optimization problem explores the possibility of “rotating” R into Q through a 
transformation matrix P which is orthonormal with respect to M. This optimization 
problem is also closely related to the following constrained algebraic problem: 

Q-RP = 0 

P T MP = /„ ( 12 ) 

We will show later that identification of linear systems with symmetric matrices can be 
formulated as solving the above constrained algebraic equations. Because || Q-RP ||2 > 0 
for all P and equality holds only when Q-RP = 0, the smallest possible minimum value 

of J is 0. Therefore, every P that satisfies Eq. (12) is also a solution of the optimization 
problem in Eq. (11). 

The derivation for a solution of the constrained optimization problem in Eq. (11) 

follows the derivation in Ref. [2] for a similar constrained optimization problem. First, J 
can be rewritten as 

J = trac e(Q T Q) + trace(P 7 ’ R T RP) - 2trac e(P T R T Q) ( 13 ) 

Note that if P satisfies the constrained condition P T MP = J n , then P t M?MiP = / 

and so, M*PP t M 2 = I n . The second term in the Eq. (13) can be converted to a form 
that is independent of P. 

tra ce(P T R T RP) = tra ce(P T M^M~^R T RM~^M^P) 

= tr<Lce(M~5 R T RM~$ M? PP T M*) 

= tra,ce(M~5 R T RM~%) 

Therefore, from Eq. (13) it is seen that minimizing J is equivalent to maximizing 

tia£e(P T R T Q). The maximizing P can be obtained by using the singular value de- 
composition technique. Let the SVD of M~tR T Q be 

M~^R t Q = UEV t 
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Since both rank(Q) and rank(-R) are < n, the rank of R T Q is also < n. Assume R T Q 
has r nonzero singular values (r < n and r is not necessarily equal to p) and arrange the 

decomposition such that 

E = diag((Ti, . . . ,<r r ,0,. . .) , &i > 0 

Define an S matrix by 

s T = v t p t m'*u 

The constrained condition P T MP = /„ implies that S is orthogonal, i.e., S T S = /„. 

Then, r 

trace! P T R T Q) = trace(P T M^f7EV T ) = trace(S T E) = ^s ti cr, 

1 = 1 

in which s„, i = 1, . . . ,p, are diagonal elements of S. Because S is an orthogonal matrix, 

su £ 1. Therefore, the maximum value of tra ce(P T R T Q) is which is attained by 

~~~ 1 = 1 

setting 

S= C with SlSb = I n -r ( 14 ) 

U Ofc 

The above derivation shows that the solution of the constrained optimization prob- 
lem in Eq. (11) can be expressed by 

P = M~'*USV t (15) 

where 5 is given by Eq. (14) and U and V are SVD of M~$R T Q. Obviously, if r < n, 
then the solution is non-unique because every S matrix in the form of Eq. (14) yields 
a minimum. Therefore, Eq. (15) represents a “solution group”, which includes all the 
possible solutions. Although every P in the solution group minimizes J, there is only 
one minimum value of J , that is 

r 

Jntin = trac e(Q T Q) + trac e{M~? R T RM*) -2^^ 

t= 1 

The solution group given in Eq. (15) is also the solution of the constrained algebraic 
problem in Eq. (12). This statement is deduced from the following arguments. As men- 
tioned previously, if there exists a P that satisfies Eq. (12), then it must be a solution of 
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the optimization problem in Eq. (11) because J =|| Q-RP ||£> 0 for all P with equality 
holds only when Q - RP = 0. All the solutions of the optimization problem are included 
by Eq. (15). If the constrained algebraic problem has any solution at all, the solution 
must also be included by Eq. (15). Since every solution in the solution group produces 
the same minimum, if one solution yields = 0, then every solution in the solution 
group yields = 0. In other words, the constrained algebraic problem in Eq. (12) 
either has no solution or has the solution given by Eq. (15). When the constrained al- 
gebraic problem has solution, the solution is unique if and only if rank(ir! r Q)=u, which 
implies rank(Q)=rank(i?)=n. In the following sections, we will use this conclusion to 
discuss identifiability of systems with symmetric matrices. 

Identification of Symmetric System Matrix 

A linear system with symmetric system matrix has real eigenvalues and a set of 
orthogonal, real eigenvectors. By using eigenvalue decomposition, any given test-data- 
denved model (A,B,C) can be transformed into the modal coordinates representation 
(A,B,C), where A € R nXn is a diagonal matrix with system eigenvalues (or poles) on 
the diagonal. Therefore, identification of A matrix in physical coordinates can be accom- 
plished from finding the transformation matrix that relates the modal-coordinates model 
and the physical-coordinates model. The modal-coordinates model and the physical- 
coordinates model are related by the following congruence transformation 

A = P T AP , B = P t B , C = CP (16) 

where P is called the modal matrix and satisfies P T P = I n . 

From the last two equations in Eq. (16), it is seen that the modal matrix P can 
be determined from solving the following constrained algebraic problem 



P T P = In 
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If P is uniquely determined, then the identification of A is accomplished by setting 
A = P.\P T . From the results derived in the previous section, P is unique if and only if 


both 
plus t 


B T 1 

C 


and 


B t 

c 


have rank equal to n. In other words, if the number of actuators 
he number of sensors is less than the number of the states, then the system matrix 


in physical coordinates is not identifiable. 


Identification of Undamped Structural Dynamics Systems 


Now, we will discuss identifiability of undamped structural dynamics systems with 
symmetric mass and damping matrices. Assume the testing environment is perfectly 
noise free. Then, the test data of an undamped structural dynamics system can be 
realized by the following modal model 

rj + = Fu 1 7 € u <E R' (18) 

y — CdV d" + C'a 1 ? V € R 

where fl 2 = diag [w?] and u>, are system’s natural frequencies. In physical coordinates, 
the system equation is described by 

Mx + Kx = Fu x € R n , ue R l 

y = Cjx + C v x + C a x y € R m 

It is assumed that physical location of actuators and sensors are given. That is, matrices 
F, C d , C v , and C a are known. We are interested in identifying M and K in physical 
coordinates from using the given modal model in Eq. (18) and the given physical locations 
of actuators and sensors. 

These modal-coordinates model and the physical-coordinates model are related by 
a modal matrix P such that 

p t mp = i n , P T KP = n 2 
F = P T F , Cd = C d P > C v = C V P , C a = C a P 

First, consider the case when the mass matrix M is given. Then, the transformation 
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matrix P can be determined from solving the following constrained algebraic problem: 


F T - 


■ F T 1 

c d 


c d 

c v 


c v 

. C a . 


Ca . 


P t MP = I n 

^ P is uni< l uel y determined, then the physical stiffness matrix is determined from setting 
K = Previous discussion shows that the total number of actuators and 

sensors must be at least equal to the number of degrees-of-freedom of the system in order 
for the modal matrix to be uniquely determined. 

For the case that both M and I< are unknown, the only equation available for 
determining P is 

■ ft - 
Cd 

C v 

C a J 

which is an unconstrained algebraic equation. After P is obtained, M and K can be 
calculated by 

m = p~ t p~ i , K = p~ T n 2 p - 1 

In order for a unique P to exist, the matrix formed by the input and output matrices 
must be consistent and must have rank n. Therefore, the knowledge about M does not 
relieve the number requirement of actuators and sensors. 

Which Part of the System is Uniquely and Exactly Identified? 

We have shown that if the number of actuators and sensors are not enough, then 
it is never possible to identify the physical system matrices uniquely. Nevertheless, it 
is still of interest to know if it’s possible to uniquely identify certain part of the system 
matrices by using the knowledge about physical locations of actuators and sensors. 

First, consider the case with unsymmetric system matrices. The physical-coordinates 
model (A, B, C) and the test-data-derived model (A, B, C) are related by a transforma- 
tion matrix P as shown in Eq. (2). The transformation matrix P is determined from 



F 1 


c d 


c v 


Ca . 


( 21 ) 
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solving Eqs. (2b) and (2c). Previous results shows that unless either rank(5)=n or 
rank(C)=n, P can not be uniquely determined. Now, assume there is only one sensor 
and it is located at the i - th physical coordinate. In this case, C = e„ where e, the t-th 
row of the identity matrix /„. Then, the relation C = CP = e t P uniquely determines 
the i-th row of P matrix. Similarly, if we assume that there is only one actuator and it is 
located at the j-th physical coordinate, then the relation B = P~ X B = P~ l eJ uniquely 
determines the j-th column of P~\ With the i-th row of P and the j-th column of P~' 
being uniquely determined, the (i, j)-th element of A matrix is also uniquely determined, 

since A = P AP " 1 . 

For systems with symmetric matrices, the test data acquired in a noise-free envi- 
ronment can be realized by a modal model (A ,B,C). The physical-coordinates model 
{A,B,C) and the modal-coordinates model (A,B,C) are related by a modal matrix P 
as shown in Eq. (16). Because P~ l = P T , either a sensor or an actuator suffices to 
determine one row of the modal matrix uniquely. A pair of noncolocated sensor and 
actuator can identify two rows of the modal matrix, and therefore, four elements of the 

system matrix. 

In summary, from knowing the physical locations of actuators and sensors, it is 
possible to identify part of the system matrix in physical coordinates. However, this 
conclusion is based on the assumption that the testing environment is noise free. In 
practical situations, the test data are inevitably contaminated by noises and small errors. 
Sensitivity of the identification results to noises is another topic worthy of studying. 

A Simple Example 

As an example we will consider the simple mass-spring system shown in Fig. 1. 
There is a displacement sensor and a force actuator located at D0F2 and D0F3, denoted 
by s and a, respectively. The modal model of this system is given by 



' 0.2708 

0 

0 


' 0.4319 ) 

V + 

0 

2.1315 

0 

» 1 = * 

-0.5524 

0 

0 

5.8477 _ 


0.0917 J 

y = 1 

[ 0.3539 

0.2325 

-0.2659 

]v 
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Assume that the mass matrix 


' 2 0 0 ' 

M = 0 4 0 

0 0 2 

- 

is given. And, we want to determine the physical stiffness matrix K. 

The modal matrix P is determined by solving the following constrained algebraic 
problem: 


' 0.4319 

-0.5524 

0.0917 ' 

0 0 1 


0.3539 

0.2325 

-0.2659 J “ | 

P 

0 10’ 

P t MP = I 3 


Or, 

Q = RP , P t MP = I 3 

Singular value decompositions of M~$R T Q(= USV) shows that there are two nonzero 
singular values, while the system has three degrees-of-freedom. Therefore, the solution 
is not unique. The solution group is given by P = M~'*USV T for all S that satisfy 

' 1 0 0 ' 

S= 0 10 , Si = \ 

.0 0 S* _ 

There are two solutions: S„ = 1 and S t = -1. Therefore, there are two modal matrix 
that are solutions of the constrained algebraic problem. They are: 

0.2510 0.2945 0.5918 1 T -0.2510 -0.2945 -0.5918 ' 

Pi = 0.3539 0.2325 -0.2659 , P 2 = 0.3539 0.2325 -0.2659 

- °- 4319 -0.5524 0.0917 J 0.4319 -0.5524 0.0917 

Using these two modal matrices to determine the stiffness matrix, we get 

9 —6 0 9 6 0 

I<i =-6 9-3 , I< 2 =6 9-3 

. 0 -3 3 0 -3 3 



Figure 1: A three degrees-of-freedom mass-spring system. 
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in which K\ is the correct stiffness matrix in physical coordinates. K 2 is not the correct 
stiffness matrix, but the (2,2), (2,3), (3,2), and (3,3) elements are correct. This is because 
that there is an actuator located at D0F3 and a sensor located at D0F2. 

This example shows that with only one actuator and one sensor (non-colocated), 
the physical stiffness matrix of the three degrees-of-freedom mass-spring system can not 
be uniquely identified from the given modal model and the given physical locations of 
actuator and sensor. However, those elements in the stiffness matrix corresponding to 
the actuator and sensor locations are exactly and uniquely determined 

Concluding Remarks 

Identifi ability of linear, time invariant systems from test data with a limited num- 
ber of actuators and sensors is investigated in this paper. It is shown that for systems 
with unsymmetric system matrices, either the number of actuators or the number of sen- 
sors must be greater than or equal to the number of states so that the system matrices 
in physical coordinates can be uniquely identified. For systems with symmetric system 
matrices, the number of actuators plus the number of sensors must be greater than or 
equal to the number of states in order for the physical system matrices to be uniquely 
determined. For undamped structural dynamics systems with symmetric mass and stiff- 
ness matrices, the total number of actuators and sensors must be greater than or equal 
to the number of degrees-of-freedom, so that the physical mass and stiffness matrices can 
be identified from test data. A simple mass-spring system is used as an example to verify 
the conclusions. 
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